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ABSTRACT 



We introduce the spatial correlation function Cq(t) and temporal autocorrelation function 
Cq{€) of the local tetrahedral order parameter Q = Q(r,t). Using computer simulations of 
the TIP5P model of water, we investigate Cq(t) in a broad region of the phase diagram. 



First we show that Cq{t) displays anticorrelation at r « 0.32nm at high temperatures 
T > TV « 250 K, which changes to positive correlation below the Widom line TV. Further 

we find that at low temperatures Cq{t) exhibits a two-step temporal decay similar to the self 

l> ■ 

intermediate scattering function, and that the corresponding correlation time tq displays a 

00 ■ 

dynamic crossover from non-Arrhenius behavior for T > TV to Arrhenius behavior for T < 

TV. Finally, we define an orientational entropy Sq associated with the local orientational 

order of water molecules, and show that tq can be extracted from Sq using an analog of the 
d ' 

Adam-Gibbs relation. 



The local structure around a water molecule arising from the vertices formed by 



est neighbors is approximately tetrahedral. The degree o 



tetrahedrality 



our near- 
[l|, bl, bi can be 



quantified by the local tetrahedral order parameter Q j4|, |5|, |y, 

Q suggest a high degree of orientational heterogeneity. The distribution of Q is bimodal, 



3, |S|] . Spatial fluctuations in 



centered around less tetrahedral and more tetrahedral values of Q. Upon decreasing tem- 
perature, the high-tetrahedrality peak grows, suggesting that the local structure of water 
becomes much more tetrahedral at low enough temperatures 9j, lid ]. 

Water has been hypothesized to belong to the class of polymorphic liquids, phase 
separating — at sufficiently low temperatures and high pressures — into two distinct liquid 
phases: a high density liquid (HDL) with smaller Q and a low density liquid (LDL) with 
larger Q [ll| . The co-existence line separating these two phases terminates at a liquid-liquid 
(LL) critical point. The locus of maximum correlation length in the one-phase region is 
called the Widom line Tw = Tw(P), near which different response functions such as iso- 

baric heat capacity Cp and isothermal compressibility Kr display maxima. Recent neutron 

I I I I 

scattering experiments [12| and computer simulations ]13] show that the dynamics of water 

changes from non-Arrhenius for T > Tw to Arrhenius for T < Tw- 
in this paper, we ask how the orientational order and its spatio-temporal correlations 
change upon crossing the Widom line. To this end, we introduce the spatial correlation 
function Cq(t) and the temporal autocorrelation function Cq(t) for the local orientational 



order o 
silicon 



water, which can also be app 



14 



ied to study other locally tetrahedral liquids such as 



15] , silica 16| , phosphorus 17| ; all of which have also been found to display some 



of the water-like anomalies. To this end, we perform molecular dynamics (MP) simulations 

9], which exhibits 
20| | . We carry out 



8. 



19 



of N = 512 waterlike molecules interacting via the TIP5P potential 

a liquid-liquid (LL) critical point at T c « 217 K and P c « 340 MPa 

simulations in the NPT ensemble at atmospheric pressure (P = 1 atm) for temperatures T 

ranging from 320 K down to 220 K. 

To quantify the local degree of orientational order, we calculate the local tetrahedral 
order parameter |4| 

s_i. _L r n 2 

(i) 



COS Ipikj + - 



where ipn-j is the angle formed by the molecule k and its nearest neighbors i, and j. The 
average value (Q) = (l/N)J2kQk increases with decreasing T, and saturates at lower T, 









while \d(Q)/dT\ has a maximum at the Widom line Tw ~ 250 K [9|, Il0| . 

To characterize the spatial correlations of the local order parameter we find all the pairs 
of molecules i and j whose oxygens are separated by distances belonging to the interval 



[r — Ar/2, r + Ar/2]. The number of such pairs is 

N(r,Ar) = ^25{r ij -r,Ar), (2) 

where r^- is the distance between the oxygens of molecules % and j. The sum is taken over 
all molecules in the system and 

( 1 if |r — TV,- 1 < Ar/2 
5(r^-r,Ar) = ' ljl ' (3) 

1 otherwise. 
For Ar — > 0, <5(rjj — r, Ar) — > 5(ry — r)Ar where 5(rjj — r) is the Dirac (^-function. iV(r, Ar) 
can be approximated as 

N(r, Ar) = Airr 2 Ng 00 {.r)p/\r. (4) 

where iV is the total number of molecules in the system, goo{ r ) is the oxygen-oxygen pair 
correlation function, and p is the number density. 

Next we find the local order parameters Qi and Qj for all such pairs of molecules and 
compute their mean 

Q(r) = (Q)r = ^V a\ - > & 

J\ [Tij — r, Ar) 

their variance 

)2 



a 2 n(r) 



(Q 2 ) r - {Q)l = £ %fi ffi - (g)?, (6) 

^rr' iV [Tij — r, Ar) 

and covariance 

inn\ = V QjQj s ( r v ~ r > Ar ) / 7 x 



'■J 



Nfrij — r ; Ar ) 



Finally we introduce the spatial correlation function of the local tetrahedrality Q 

Co(r) = — w> — ' ( ' 

In Figs.[T](a) andQJb), we show Q(r) and (y%{r) for different temperatures and atmospheric 
pressure. The behavior of Q(r) and its variance ^{r) as a function of the distance r has a 
clear physical meaning. For the molecules separated by 0.32 nm, Q(r) has a deep minimum 
characterizing the distortion of the first tetrahedral coordination shell of the four nearest 
neighbors by the intrusion of a "fifth neighbor" 21j. Conversely, the quantity <Jn{r = 



0.32nm) — Oq(oo) for such molecules dramatically increases upon decreasing temperature [see 
Fig. Q](c)], and has a maximum at T « 246 K which is approximately equal to the ambiant 
pressure value of T w (« 250 K for the TIP5P model) reported in Refs. [lOj, [l3|. In Fig. [D(d), 



we show Cq(t) for different temperatures for P = 1 atm. Cq(t) has positive maxima 
at positions of the first and second peaks in the oxygen-oxygen pair correlation function, 
suggesting that the molecules located at the tetrahedral positions are strongly correlated in 
their Q values. Water molecules separated by r « 0.32nm exhibit weak anticorrelation in 
local tetrahedral order at high temperatures, which changes to positive correlations upon 
decreasing temperature below Tw- 

To study the time development of local tetrahedral orientational order parameter, we 
introduce the time autocorrelation function 

° Q{t) = WTW ' {) 

In Fig. [2](a) we show Cg(t) for different temperatures. The decay of Cq{t) is reminiscent of 
the decay of the self intermediate scattering function. The long time behavior of Cq(£) is 
exponential at high T, but at low T can be fit with a stretched exponential exp[— (t/r)^], 
where < /3 < 1. We define the correlation time tq as the time required for Cq(t) to decay 
by a factor e. Figure 12(b) shows the values of tq as function of 1/T on an Arrhenius plot. 
The behavior of tq is non-Arrehenius at high temperatures and can be fit by a power law 



(T — Tmct) -7 where Tmct ~ 235 is the mode coupling temperature [13|, |22J. At low T, tq 
deviates from the power-law fit and becomes Arrhenius. This crossover in relaxation of local 
orientational order occurs near Tw- 

Since Q measures the local orientational order, it must contribute to the entropy of the 
system. We next derive an expression for this "orientational entropy" S(Qi, Q 2 , ■ ■-, Qn), 
which we define to be the logarithm of the number of states corresponding to the in- 
terval between (Q 1 , Q 2 , ...., Qn) and (Qi + AQi,Q 2 + AQ 2 ,-..,Qn + AQ N ). According 
to Eq. ([!]), 8(1 — Qfc)/3 = const defines a surface of a six-dimensional hypersphere of 



radius ^8(1 — Qk)/3 in the space defined by the six tetrahedral angles ipikj of Eq. (CQ). 
Hence we assume that the number of states between Qk and Qk + AQk of molecule k 
scales as (1 — Qk) AQk- We assume that the order parameters of each molecule is in- 
dependent, an assumption justified by the small value of Cg(r) studied. Thus we can 
define the number of states fl(Qi,Q 2 , ■■■■,Qn) m the interval between (Qi,Q 2 , ■■■■Qn) and 
(Qi + A<3i, Q 2 + AQ 2 , ..., Q N + AQ N ) as the product 

N 

fi(Qi, Q2, Qs, :.;Qn) = nq J|(i -Q k )K (10) 



fc=i 



where Qq = const. Hence the orientational entropy of the entire system is given by 

N N 

S(Qi, Q 2 , Q 3 , -, Qn) = NS + In JJ(1 - Q k )l = ns + -J2 M 1 - <?*)• (n) 



2 
k=i fc=i 



where So = lnO . If P(Q,T) is the distribution of Q at a given temperature T, then the 
orientational entropy Sq(T) per particle at temperature T can be written as 

S Q (T) = S + - ln(l -Q)P(Q,T)dQ. (12) 

•^ Wmin 

In Fig. [3] (a), we show that Sq(T) — So decreases with decreasing temperature as expected. 

We further define a measure of "orientational specific heat" as Cp(T) = T(dSQ(T)/dT)p. 
Figure Mjo) shows the temperature dependence of Cp(T) for P = 1 atm. Cp(T) shows a 
maximum at T « 245K rj T^y, where the total specific heat Cp has a maximum J10|, Il3| . 
suggesting that the fluctuations in orientational order reach a maximum at T\y{P)- 

To relate the orientational entropy Sq(T) to the orientational relaxation time tq(T) 
associated with the local tetrahedral ordering, we propose the following generalization of the 
Adam-Gibbs relation between the translational relaxation time and configurational entropy 

2d, y, 



r Q (T)=r Q (0)ex V [A/TS Q (T)], (13) 

where tq(0) is the orientational relaxation time at very large T, and A is a parameter playing 
the role of activation energy. Accordingly we calculate [Fig. |2](b)] Tq(T) from Sq(T) using 
Eq. ffl3|) with three free parameters: A, r(0), and Sq. We find that So ~ 4.33. 

Figure EJ^b) is an Arrhenius plot of tq, calculated from Eq. (TT3|) . The temperature 
dependence of tq is different at low and high temperatures, changing from non-Arrhenius 
(a T-dependent slope on the Arrhenius plot) which can be fit by a power law at high 



12 



Q. 



temperatures to Arrhenius (a constant slope) at low temperatures [101 , 

In summary, we have studied the space and time correlations of local tetrahedral order, 
presumably related to local orientational heterogeneities. We find that the spatial correlation 
of the the local tetrahedral order is anticorrelated for the molecules separated by 3.2A at 
high temperatures. This negative correlation changes to positive correlation upon decreasing 
T below the Widom temperature Tw{P)- Further, we define a measure of the orientational 
entropy Sq and find that it well describes orientational relaxation using the Adam-Gibbs 
relation. 
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FIG. 1: (color online) (a) The average order parameter Q(r) and (b) its variance, crg(r), as a 
function of the distance r. (c) Temperature dependence of o%{r = 0.32nm) — a% (oo) shows a 
maximum at the Widom temperature, suggesting the local fluctuations in Q at the fifth-neighbor 
distance increases upon decreasing temperature and has a maximum at the Widom line, (c) Spatial 
correlation function Cn(r) of orientational order parameter Q [Eq. ([!])], at various temperatures 
for pressures P = 1 atm. Cq{r) has positive peaks at the positions of the nearest neighbor peaks in 
oxygen-oxygen pair correlation function goo{r). A negative minimum at the fifth neighbor distance 
r ~ 0.32 nm for high T implies that the local tetrahedral order of a central molecule and its fifth 
neighbor are anticorrelated at T > 250K. Interestingly the anticorrelation at r = 0.32 nm changes 
to positive correlations upon crossing below the Widom line around T\y ~ 250 K. 
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FIG. 2: (color online) (a) Autocorrelation function Cn(t) of orientational order parameter Q at 
various temperatures. Cq(1;) is exponential at high temperatures but displays a visible two step 
decay at low temperatures, (b) Correlation time tq extracted from Cq{€). Solid line is the fit using 
the Adam-Gibbs relation [Eq. (|13|) ] between the orientational entropy Sq(T), and the orientational 
relaxation time tq. The dotted lines in (b) is the power law fit (T — Tmct)~ 7 - Behavior of tq 
deviates from the power law fit and has a crossover to Arrhenius behavior at low temperature 
near the same temperature T\y where the response functions have maxima, which we identify with 
crossing the Widom line T\y. 
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FIG. 3: (color online) Temperature dependence of (a) orientational entropy Sq(T), defined in 
Eq. (jlip . and (b) orientational specific heat Cp = T (dSQ/dT) p , which has a maximum around 
the same temperature, T ~ 246 K ~ Tyy, where the total specific heat Cp has a maximum. Solid 
line in (a) is a 5 i/l -order polynomial fit through the data which is used to generate (b) 



